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The F-GAMMA program is a coordinated efTort to investigate the physics of Active Galactic Nuclei 
(AGNs) via multi-frequency monitoring of Fermi blazars. The current study is concerned with the 
broad-band radio spectra composed of measurement at ten frequencies between 2.64 and 142 GHz. It 
is shown that any of the 78 sources studied can be classified in terms of their variability characteristics 
in merely 5 types of variability. The first four types are dominated by spectral evolution and can 
be reproduced by a simple two-component system made of the quiescent spectrum of a large scale 
jet populated with a flaring event evolving according to The last type is characterized by an 
achromatic change of the broad-band spectrum which must be attributed to a completely different 
mechanism. Here are presented, the classification, the assumed physical system and the results of 
simulations that have been conducted. 



I. INTRODUCTION 



Among the most evident characteristics of blazars 
is the intense variability at all wavelengths. Studies of 
the variability characteristics, preferably with simul- 
taneous data, can shed light on the physics driving 
the energy production and dissipation in these sys- 
tems 

[e.g-MH- The F-GAMMA program [see|j-l6| is 
a coordinated effort to explore exactly this very pos- 
sibility by monthly monitoring of Fermi blazars. F- 
GAMMA is covering mostly the radio cm to sub-mm 
bands primarily with the Effelsberg 100-m, the IRAM 
30-m and the APEX 12-m telescopes (although opti- 
cal telescopes are participating as well, Fuhrmann et 
al. in prep.) for roughly 60 prominent blazars. 

The cause for the variability itself has been long de- 
bated. The "shock-in- Jet" model suggested by is 
the most accepted one and attributes the variability 
to shocks propagating down the jet. The basic as- 
sumption is that changes at the onset of the jet, (e.g. 
changes in the injection rate, the magnetic field, bulk 
Lorentz factor etc.) cause the formation of shocks, 
which then suffer first Compton, then synchrotron and 
finally Adiabatic losses. This is the main model pre- 
diction which is used in the following. 

Here it is argued that (a) the variability patterns 
can be classified in only 5 phenomenological types and 
(b) the different phenomenological classes can be re- 
produced with a simple system composed of a quies- 
cent spectrum populated by a flaring event evolving 
according to the "shock-in- Jet" model. 



II. OBSERVATIONS AND DATA 
REDUCTION 

The observations discussed here have been con- 
ducted quasi-simultaneously with the Effelsberg 100- 
m and the IRAM 30-m telescope (the combined spec- 
tra coherency time is a few days) within the F- 
GAMMA program The 100-m telescope has 

been observing between 2.64 and 43.05 GHz at 8 fre- 
quencies and the 30-m telescope at 86 and 142 GHz 
(details in Fuhrmann et al. in prep., Nestoras et al. 
in prep., Angelakis et al. in prep.). In the current 
study only data collected until June 2011, have been 
used. 

The data reduction includes: (a) Pointing cor- 
rection, (b) Elevation dependent gain correction, (c) 
Atmospheric opacity correction, (d) Absolute calibra- 
tion (sensitivity correction). The overall uncertainties 
reached for the F-GAMMA program are of the order 
of 0.5 - 5 % for Effelsberg and of the order of < 10 % 
for IRAM. More details can be found in ]7| as well as 
in Angelakis et al. (in prep.), Fuhrmann et al. (in 
prep.) and Nestoras et al. (in prep.). 



III. ANALYSIS 

A. Phenomenological Classification of the 
Variability Patterns 

The visual inspection of the examined sources re- 
veals a plurality in spectral features as well as in the 
variability pattern that different sources exhibit. De- 
spite the apparent complexity it appears that any of 
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(a)Type 1 



(b)Type lb 



(c)Type 2 





(d)Type 3 



(e)Type 3b 



(f)Type 4 






(h)Type 5 

FIG. 1: The prototypes of the 5 types and their sub-types. 



(i)Type 5b 



the 78 sources studied here, can be classified in one 
among only five phenomenological classes on the ba- 
sis of its variability pattern, which are termed numer- 
ically from 1 to 5 (more details will be given by Ange- 
lakisetal. in prep.). Four of them show also sub- types 
which however do not deserve a separate type and are 
named after the main type followed by the letter "b" 



The prototype sources are shown in Figures 1(a) l(i) 
Their phenomenological characteristics, are: 
Type 1: is clearly dominated by spectral evolution. 
At an instant in time the spectrum appears convex 
and its peak is drifting within the observing band- 
pass from high towards lower frequencies, covering a 
significant area in the S — v space. The convex compo- 
nent is smoothly changing towards an ultimate flat or 
mildly steep power-law which is then followed by con- 
sequent events. There is no evidence for a stable steep 
spectrum. The lowest frequencies in the bandpass are 
remarkably variable indicating that the activity seizes 
at frequencies much lower than the lowest in our band- 



pass. The prototype source is shown in Figure 1(a) 
Type lb: As a sub-class of the previous one, type lb 
shows similar characteristics except that the lowest 
frequency does not show as intense variability. The 
activity seizes around this part of the band-pass (see 
Figure 1(b) ). 



Type 2: is also dominated by spectral evolution. The 
basic characteristic of this case is the fact that the 
flux density at the lowest frequency during the steep- 
est spectrum phase is higher than that during the in- 
verted spectrum phase. Moreover, the maximum flux 



density reached by the flaring events is significantly 
above that at the lowest frequency. This implies that 
the observed steep spectrum is not a quiescent spec- 
trum but rather the "echo" of an older, yet recent, 
outburst. The prototype of this type is shown in Fig- 
ure 1(c) 

Type 3: Type 3, shown in Figure l(d)[ is also dom- 
inated by spectral evolution. The identifying charac- 
teristics of this type are: (a) the fact that the lowest 
frequency practically does not vary and, (b) the maxi- 
mum flux density level reached by outbursts is compa- 
rable to that at the lowest band-pass frequency. This 
phenomenology leaves hints that the events seize very 
close to the lowest frequency of the band-pass and 
hence a quiescent spectrum is becoming barely evi- 
dent. 

Type 3b: Type 3b, shown in Figure 1(e) is very sim- 
ilar to type 3. Here however the quiescent spectrum 
is seen clearly at least at the 2 lowest frequencies. 
Type 4: Sources of this type spend most of the time 
as steep spectrum ones which are sometimes showing 
an outburst of relatively low power propagating to- 
wards low frequencies. A representative case is shown 



in Figure 1(f) 

Type 4b: This type includes persistently steep spec- 



trum cases as it is shown in Figure 1(g) 



All previous classes, are clearly dominated by spec- 
tral evolution. There exists a class of sources for which 
the variability happens self-similarly without signs of 
spectral evolution. Those are grouped in a separate 
type with two sub-types: 
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Type 5: In this case the spectrum is convex and fol- 
lows an "achromatic" evolution. That is, it shifts its 
position in the S — v space preserving its shape. This 



is shown clearly in Figure 1(h) 



Type 5b: This type shows, in principle, character- 
istics similar to the previous one but there occurs a 
mild yet noticeable shift of the peak (5'm, i^m) towards 
lower frequencies as the peak flux density increases. 



A characteristic case is shown in Figure l(i) 



This classification is done solely on the basis of 
the phenomenological characteristics of the variabil- 
ity pattern shown by the radio spectra within a given 
band-pass. As it is discussed in the next section, it 
appears that all the phenomenology for types l-4b 
can be naturally explained with the same underlying 
system observed under different circumstances. 



B. A Physical Interpretation of the Variability 

Types 1 - 4b 

The phenomenological types 1 - 4b discussed earlier 
can be reproduced by the same simple two-component 
system, made of: (a) A power-law quiescent spec- 
trum attributed to the optically thin emission (large- 
scale jet and recent flaring events) and (b) the convex 
synchrotron self-absorbed spectrum of a current out- 
burst superimposed on the quiescent part. The as- 
sumed configuration is presented in Figure [2] where 
the shaded areas denote the observing band-pass. 
The phenomenology shown there captures the system 
(solid line) at an instant in time and the spectral shape 
that would be observed depends on: (a) the position 
of the shaded areas relative to the high and low fre- 
quency peak (i.e. the peak of the outburst) and (b) 
the width of the band-pass relative to the width of 
the bridge between the optically thick part of the out- 
burst and the steep part of the quiescent spectrum. 

These two quantities can be modulated by the com- 
bination of (a) the redshift and (b) the source in- 
trinsic properties. The redshift changes the rela- 
tive position of the band-pass allowing a different part 
of the spectrum to be sampled. The source intrinsic 
properties imply that different sources show different 
spectral characteristics (e.g. peak frequency of the 
outburst, peak flux density excess of the outburst over 
the quiescent spectrum, broadness of the valley etc.). 
More importantly, the dynamical evolution of a flaring 
event, is a function of the source intrinsic properties 
and introduces a third parameter (c) the flare specific 
properties which is determining the characteristics of 
the variability pattern. 

In order to examine whether the assumed model 
can reproduce the observed phenomenologies, some 
characteristic cases have been evaluated. 




Gl 



FIG. 2: The assumed two-component sys- 
tem. The different variabihty types can 
be reproduced with the appropriate mod- 
ulation of the relative position and relative 
broadness of the band-pass denoted by the 
grey shaded areas. 



Reproducing the Observed Phenomenologies 
1 - 4b 



Following the hypothesis that all the observed 
events are the reflection of the same process, namely 
shocks evolving in jets seen with different frequency 
band-passes at different evolutionary stages, the 
shock-in-jet model [l,, ,8ij has been applied to reproduce 
their temporal evolution. The followed approach is 
presented in Q where the flaring event passes different 
radiative evolutionary stages {Compton, Synchrotron 
and Adiabatic stage). Simulations have been made for 
a large fraction of the parameter space. In Figure [3] 
are presented only three cases of sources at z — 1.5 
but of different luminosity. From these plots it is al- 
ready clear that the assumed scenario can reproduce 
most of the observed phenomenologies. 



IV. DISCUSSION 

The variability patterns of the studied blazars can 
be categorized in (a) spectral evolution dominated 
cases and (b) self-similarly varying convex spectra 
ones. This implies that there must exist two distinct 
mechanisms causing variability. This refers to the 
available baseline (roughly five years) meaning that 
sources of type 5 and 5b could still show spectral evo- 
lution over longer time scales. An additional element 
that points towards a different variability mechanism 
is the persistency of the spectral shape, in the case of 
type 5 and the fashion of change in the case of type 5b 
which is not seen in cases of clear spectral evolution. 

Of the 78 examined sources, 8 show achromatic 
variability. The interesting characteristics is that the 
cases that show a mild spectral evolution (type "5b"), 
the turnover flux and frequency Sm and t'm , are evolv- 
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FIG. 3: Three characteristic cases of the system assumed for the reproduction of the phenomenologies of types l-4b. 



ing in an anti-correlated fashion (see e.g. Figure l(i) ) 



Apart from the fact that all of them show the clear 
presence of a large scale jet even at 2 cm (as it is shown 
form the MOJAVE images Kellermann et al. 10) no 
other peculiar property has been identified so far. Pos- 
sible mechanisms that are examined to be producing 
this variability include: opacity effects, changes in the 
magnetic field structure, changes in the Doppler fac- 
tors and geometrical effects. 

None of the studied sources has shown a switch of 
type over the baseline of the F-GAMMA program, nei- 
ther between types of the same underlying mechanism 
(i.e. l-4b) nor between types with different underly- 
ing mechanism (i.e. types l-4b and types 5, 5b). This 
suggests that the mechanisms producing the variabil- 
ity is either a fingerprint of the source or the condi- 
tions that determine it change over longer time scales. 
It would be essential to observe whether a source can 
switch from achromatic to an evolution dominate type 
or even further whether it exhibits periods of either 
variability flavor. In any case, the persistency of the 
evolution dominated types implies that the power de- 
posited in each event for a certain source is not vary- 
ing significantly from one event to the other. Further 
investigations to investigate this statement are under- 
way and will be presented elsewhere. 

Concerning the evolution dominated case, it seems 



that the model provides a precise reproduction of 
the observed phenomenology and most importantly, 
over a range of intrinsic parameters covered by the F- 
GAMMA sample. Studies to examine whether other 
variability mechanisms can reproduce the observed 
phenomenology (shapes, time scales etc.) are needed. 
In any case, any successful model could be used for 
reversing the process and be used for calculating phys- 
ical parameters from the observed spectra. 
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